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The geometric frustration of the spin-1/2 Ising-Heisenberg model on the triangulated Kagome 
(triangles-in-triangles) lattice is investigated within the framework of an exact analytical method 
based on the generalized star-triangle mapping transformation. Ground-state and finite-temperature 
phase diagrams are obtained along with other exact results for the partition function, Helmholtz free 
energy, internal energy, entropy, and specific heat, by establishing a precise mapping relationship to 
the corresponding spin-1/2 Ising model on the Kagome lattice. It is shown that the residual entropy 
of the disordered spin liquid phase is for the quantum Ising-Heisenberg model significantly lower 
than for its semi- classical Ising limit {So/NTkB ~ 0.2806 and 0.4752, respectively), which implies 
that quantum fluctuations partially lift a macroscopic degeneracy of the ground-state manifold in the 
frustrated regime. The investigated model system has an obvious relevance to a series of polymeric 
coordination compounds Cu9X2(cpa)6 (X=F, CI, Br and cpa=carboxypentonic acid) for which we 
made a theoretical prediction about the temperature dependence of zero-field specific heat. 
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I. INTRODUCTION 



The antiferromagnetic quantum Heisenberg model 
(AF-QHM) defined on geometrically frustrated planar 
lattices represents a long-standing theoretical chaUenge 
due to a rich variety of unusual ground states it exhibits 
as a result of the mutual interplay between quantum fluc- 
tuations and geometric frustrationii^i^i^ In particular, 
the extensive theoretical studies of the spin-1/2 AF-QHM 
on the triangular;^ Kagome^^ Shastry-Sutherland^ 
star,^ checkerboard^^ square Kagome latticeii and 
others have revealed a great diversity in their ground- 
state and low-temperature behavior. It is now widely 
accepted that the ground state of the spin-1/2 AF-QHM 
on some geometrically frustrated planar lattices like a 
triangular lattice is the Neel-like ordered stated while 
there is a still controversial debate whether— or noii the 
disordered spin liquid state is the true ground state of 
this model on Kagome lattice. Anyway, the macroscopic 
degeneracy of the ground-state manifold turns out to be 
strongly related to a geometric topology of the underlying 
lattice and it is therefore of particular research interest 
to explore a connection between the zero-point entropy 
and the lattice geometry. 

Another striking feature, which currently attracts a 
great deal of attention to the spin-1/2 AF-QHM on the 
geometrically frustrated planar lattices, is being a pres- 
ence of quantized magnetization plateaux in their low- 
temperature magnetization curves. ^''^ It is worthwhile to 
remark that this outstanding quantum phenomenon has 
already been experimentally observed in several proto- 



typical examples of the frustrated quantum spin systems 
such as the triangular lattice compounds CsFe(S04)2r^ 
Cs2CuBr4fi^ and RbFe(Mo04)2r^'^ the Kagome lattice 
compound [Cu3(titmb)2(CH3COO)6].H20,i6 and the 
Shastry-Sutherland lattice compounds SrCu2(B03)2 — 
and RB4 (R = Er, Tm).^* It also has been demonstrated 
that the interplay between the geometric frustration and 
quantum fluctuations might be a driving force for an en- 
hanced magnetocaloric effect emerging during the adia- 
batic demagnetizationji^ This makes from geometrically 
frustrated spin systems especially promising refrigerant 
materials in view of reaching temperatures in a sub- 
Kelvin range, since they often remain disordered down 
to the lowest achievable temperatures unlike paramag- 
netic salts usually exhibiting a spin-glass transition. 

The aforementioned scientific achievements have stim- 
ulated an intensive search for inorganic molecular ma- 
terials, whose paramagnetic metal centres connected in 
the crystal lattice via superexchange pathways would be 
strongly frustrated by their geometric arrangement 3^ 
From this perspective, the series of isostructural 
polymeric coordination compounds Cu9X2(cpa)6.nH2 
(X=F, CI, Br and cpa=carboxypentonic acid)£L belongs 
to the most fascinating geometrically frustrated materi- 
als on behalf of a beautiful architecture of their magnetic 
lattice. The magnetic lattice of this series is built up of di- 
valent copper (Cu^+) ions situated at two crystalograph- 
ically inequivalent lattice positions (see Fig. [1]). Cu^+ 
ions with a square pyramidal coordination (a-sites) form 
equilateral triangles (trimers), which are inter-connected 
to one another by Cu^^ ions (monomers) with an elon- 
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gated octahedral coordination (6-sites). The lattice po- 
sitions of the &-sites constitute the regular Kagome net- 
work, whereas each monomeric &-site is connected via 
four bonds to two adjacent trimers of the a-sites. This 
magnetic architecture can be accordingly regarded as 
the triangulated Kagome (triangles-in-triangles) lattice, 
since smaller triangles of the trimeric a-sites are in fact 
embedded in larger triangles of the monomeric 6-sites 
forming the basic Kagome pattern. 

Experimental studies reported on this family of com- 
pounds reveal obvious manifestations of the geometric 
frustration. All three isomorphous compounds do not or- 
der down to 1.7 the magnetization shows a plateau 
around one third of the saturation magnetization and it 
does not saturate up to 38 T.— Furthermore, the temper- 
ature dependence of inverse susceptibility indicates two 
temperature regimes inherent to two different exchange 
pathways. The linear dependence of the inverse suscep- 
tibility is well fitted by the Curie- Weiss law within the 
temperature interval between 150-250 K with the respec- 
tive Weiss constants = -237 K, -226 K, and -243 K 
for the fluoro, chloro, and bromo analogue^, respectively, 
while the Weiss constant of the best Curie- Weiss fit gen- 
erally increases up to roughly = 6 K in the tempera- 
ture range below 50 K.^"^ These observations would sug- 
gest an extremely high frustration ratio / — \Qv,\/Tc > 
130 (Tc is the ordering temperature)^ for each mem- 
ber of the CugX2(cpa)6 family, which makes from this 
series a prominent class of the highly frustrated mate- 
rials that possibly display the spin liquid ground state. 
Based on the considerations of exchange pathways, the 
strong antiferromagnetic interaction has been assigned 
to the exchange interaction Jaa between the trimeric a- 
sites, while the weaker (possibly ferromagnetic) interac- 
tion has been ascribed to the exchange interaction Jat 
between the trimeric a-sites and the monomeric 6-sites. 
The overall ratio between both exchange constants might 
be estimated from the corresponding ratio between the 
Weiss constants extrapolated from the high- and low- 
temperature regime yielding \Jaa/Jab\ ~ 40. 

Motivated by these experiments, Zheng and Sun^^ 
have calculated an exact phase diagram of the spin-1/2 
Ising model on the triangulated Kagome lattice and the 
validity of this phase diagram has recently been con- 
firmed by an independent exact calculation of Loh, Yao, 
and Carlsoni^ In addition to the exact ground-state and 
finite-temperature phase diagrams, the authors of the lat- 
ter work also presented exact analytical results for several 
thermodynamic quantities (partition function, Helmholtz 
free energy, internal energy, entropy, and specific heat), 
which were complemented by the Monte Carlo simula- 
tions corroborating these exact analytical results and 
bringing other accurate numerical results for the mag- 
netization and susceptibility in the zero as well as non- 
zero magnetic field.-^S Among the most interesting find- 
ings being reported is certainly a theoretical prediction 
of the spin liquid phase with a large residual entropy per 
spin ^o/A^T^B = ^lii'''2 = 0.4752, which appears in the 



ground state on assumption of a sufficiently strong an- 
tiferromagnetic intra-trimer interaction Jaa/\Jab\ < — 1 
(see for details Section IV. A in Ref. [26l ). 

Unfortunately, the theoretical description based on the 
Ising model might fail in describing many important 
vestiges of the most of copper-based coordination com- 
pounds, since they usually exhibit a rather isotropic mag- 
netic behavior. Beside this, the Ising model description 
completely disregards the effect of quantum fluctuations, 
which might play a crucial role in determining the mag- 
netic behavior of coordination compounds incorporat- 
ing paramagnetic Cu'^^ ions having the lowest possible 
quantum spin number S = 1/2. Therefore, the spin-1/2 
quantum Hcisenbcrg model is usually thought of as be- 
ing much more appropriate model for the copper-based 
coordination compounds In accordance with this state- 
ment, ESR measurements performed on the Cu9X2(cpa)6 
series have shown just a minor anisotropy in the (7-factor 
that also serves in evidence of a negligible magnetic 
anisotropy 

Of course, the spin-1/2 quantum Heisenberg model 
is very difficult to deal with due to insurmount- 
able mathematical complexities associated with a non- 
commutability between spin operators involved in its 
Hamiltonian and thus, one has usually to rely either 
on applicability of some simpler approximative method 
or to perform rather extensive numerical calculations. 
To the best of our knowledge, the spin-1/2 Heisen- 
berg model on the triangulated Kagome lattice has 
been studied yet merely in terms of the linear Holstein- 
Primakoff spin wave theory^^ and the variational mean- 
field like treatment.'^" Both the aforementioned methods 
have however obvious insufficiencies. The former method 
based on the spin- wave approximation is applicable only 
if the ratio between both interaction parameters is from 
the interval 1 < \ Jaa/Jab\ ^ 3, i.e. the frustrating anti- 
ferromagnetic intra-trimer interaction Jaa might not be 
much too stronger than the ferromagnetic one Jab- On 
the other hand, the latter method based on the varia- 
tional mean-field like treatment is asymptotically exact 
in the \Jab/Jaa\ — ^ limit, but it even fails to predict the 
disordered ground state for any finite ratio \Jab/ Jaa\^ 

In view of this, the present work aims to suggest and 
exactly solve the spin-1/2 Ising-Heisenberg model on the 
triangulated Kagome lattice by establishing a precise 
mapping relationship with the corresponding spin-1/2 
Ising model on the simple Kagome lattice. Within the 
framework of the proposed Ising-Heisenberg model, the 
exchange interaction between the trimeric a-sites is being 
treated as the XXZ Heisenberg interaction Jaa = •/h(A), 
while the exchange interaction between the trimeric a- 
sites and monomeric 6-sites will be approximated by the 
Ising- type interaction Jab = Ji- Even although this 
model also has a clear deficiency in that the monomer- 
trimer interaction is being considered as the Ising-type 
interaction, it should be much more reliable in describ- 
ing a frustrated magnetism of CugX2(cpa)6 compounds 
as it correctly takes into account quantum fluctuations 
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FIG. 1: The cross-section of the triangulated Kagome (triangles-in-triangles) lattice. The empty and full circles denote lattice 
positions of the trimeric o-sites and monomeric 6-sites, respectively, which are occupied within the proposed Ising-Heisenberg 
model by the Heisenberg and Ising spins. Solid and broken lines schematically reproduce the intra-trimer Heisenberg interaction 
Jaa = Jii{^) and the monomer-trimer Ising interaction Jab = Ji- The ellipse demarcates a six-spin cluster described by the 
Hamiltonian ([2| , which can be mapped by the use of generalized star-triangle transformation into a simple triangle of the Ising 
spins mutually interacting via the effective exchange coupling /3J^fg. 



between the trimeric a-sites. In addition, it is also quite 
plausible to argue that the monomeric &-sites are in the 
frustrated regime completely free to flip without any en- 
ergy cost and hence, the Ising character of the monomer- 
trimer interaction might be regarded as at least quite 
reasonable first-order approximation. 

The rest of this paper is organized as follows. In Sec- 
tionHH we will provide a detailed description of the model 
under investigation and then, basic steps of our exact an- 
alytical treatment will be explained. The most interest- 
ing results are presented and detailed discussed in Section 
mil where a particular emphasis is laid on a physical un- 
derstanding of the ground-state and finite-temperature 
phase diagrams, as well as, the temperature dependences 
of several thermodynamical quantities. Some conclusions 
are finally drawn in Section [IVI 



II. MODEL AND ITS EXACT SOLUTION 

Let us define the spin-1/2 Ising-Heisenberg model 
on the triangulated Kagome lattice (Fig. [1]), which 
resembles a rather curious magnetic structure discov- 
ered in the series of polymeric coordination compounds 
Cu9X2(cpa)6'^ The magnetic properties of this series are 
captured to the lattice positions of Cu^+ ions, which are 
situated at two inequivalent lattice positions previously 
referred to as the trimeric a-sites (empty circles) and the 
monomeric 6-sites (full circles), respectively. Let us now 
assign the Heisenberg spin S ~ 1/2 to each trimeric a-site 



and the Ising spin fj, = 1/2 to each monomeric &-sitej^ 
In this way, the exchange interaction Jaa between Cu^"*" 
ions located at the nearest-neighbor trimeric sites will be 
treated as the Heisenberg interaction Jh(A), while the 
exchange interaction Jab between two Cu^+ ions located 
at the nearest-neighbor trimeric and monomeric sites will 
be treated as the Ising interaction Jj. The total Hamil- 
tonian of the model under investigation then reads 



n 



[A(5f^J 



SfS' 



Ji^4^Af,(i) 
(fc,i) 



where the first summation is carried out over all pairs 
of the nearest-neighbor Heisenberg spins and the second 
summation extends over all pairs of the nearest-neighbor 
Heisenberg and Ising spins, respectively. The spin op- 
erators Sf(a = x,y,z) and /i^ denote spatial compo- 
nents the usual spin-1/2 operator and the parameter A 
will allow us to control the exchange anisotropy in the 
anisotropic XXZ Heisenberg interaction and to obtain 
the Ising model as a special limiting case for A = 0. Fi- 
nally, the total number of the Ising spins (6-sites) is set 
to N , which implies that the total number of all spins 
(lattice sites) is = 3N. 

For further convenience, the total Hamiltonian ^ can 
be written as a sum over six-spin cluster Hamiltonians 
H. = J2k'^k, whereas each cluster Hamiltonian Hk in- 
volves all the interaction terms associated with three 
Heisenberg spins from fcth trimer (see Fig. [1]) 



- Jh 



^i^kl^k2 + ^kl^k2) + Skl^k2 



- Jh 



A{S^2^k3 + ^k2^k3) + ^k2^, 



^k3 



kl) 



— S^iihki + hk2) ~ Sl2{hk2 + hks) — Sl^{hk3 + hki), 



(2) 
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where hki = ^lAL- Obviously, the six-spin cluster Hamil- 
tonian ^ might be regarded as the Hamiltonian of kth 
Heisenberg trimer placed in some generally non-uniform 
magnetic field produced by three surrounding Ising spins. 
By taking into account the commutability between differ- 
ent cluster Hamiltonians [HijHj] ~ valid for each i ^ j, 
the partition function of the spin-1/2 Ising-Heisenberg 
model on the triangulated Kagome lattice can be par- 
tially factorized into the product of cluster partition func- 
tions Zk 

2N/3 2N/3 

^ = E n TrfcCxp(-/37T!fe) = U (3) 

{m.} k=l {m,} k=l 

Above, the summation } runs over all possible spin 



configurations of the Ising spins and the symbol Tr^ 
stands for a trace over spin degrees of freedom of kth 
Heisenberg trimer. To proceed further with the calcula- 
tion, it is very convenient to accomplish an exact analyt- 
ical diagonalization of the cluster Hamiltonian ([2]) in a 
particular Hilbert subspace corresponding to fcth Heisen- 
berg trimer. If doing so, the resulting cluster partition 
function Zk will depend just upon three Ising spins /ifei, 
/i/c2, and /^fe3, which are included in the parameters hki, 
hk2, and hks determining the 'local magnetic fields' act- 
ing on the Heisenberg spins. Moreover, the explicit math- 
ematical form of the cluster partition function Zk imme- 
diately implies a possibility of performing the generalized 
star-triangle mapping transformatioit^li^^ 



Zk = 



2exp cosh [f3Wh + f^l2 + mIs)] 
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A exp [/3J^fg (/ifeiMfe2 + Mfc2Mfc3 + MlsMfei)] : 



(4) 



where /3 = l/(fcBT), fee is Boltzmann's constant, T the 
absolute temperature, and the parameters P, Q±, and 
^± are defined as follows 



^ ~ (MfclMfc2 + Mfc2Mfc3 + MfesMfel) 



1 (Jl 



±- , 

2 V 3 



■ arctan 



'P^ - Ql 



Q± 



(5) 

;] 
(6) 

(7) 



From here onward, our procedure closely follows the ap- 
proach developed by Zheng and Sun that relates an ex- 



act solution of the spin-1/2 Ising model on the triangu- 
lated Kagome lattice to an exact solution of the the cor- 
responding spin-1/2 Ising model on the simple Kagome 
lattice.^- As a matter of fact, the star-triangle transfor- 
mation ([4]) effectively removes all the interaction terms 
associated with fcth Heisenberg trimer and substitutes 
them by the effective interaction between three en- 
closing Ising spins ^ki, Mfc2, and /i^a. Of course, the 
mapping transformation ([4]) must hold for any available 
spin configuration of three enclosing Ising spins and this 
self-consistency condition unambiguously determines so 
far not specified mapping parameters A and JJa^ 

A ^ {V^Vif and /3 jfg = In (^1/1/2) . (8) 

The functions Vi and V2 entering into the effective map- 
ping parameters A and PJ^^ are actually two indepen- 
dent expressions for the cluster partition function @ to 
be obtained by considering all eight possible spin config- 
urations available to three enclosing Ising spins 
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Vi = 2exp ( ^/MHj cosh + 2exp ( — ) cosh 
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^2 = 2 exp ( I cosh I — 



+ exp 



/3Jh /3Ji 



exp 



exp 



exp (/3JhA) + 2cxp 
-2/3sgn((7+)y^cos 



, , , 27m 

-2^sgn(g_)7^cos ( 0_ + — 



/3JhA 



27rn\ 



(9) 



(10) 



r 



and the parameters p, g± , and (fi± are defined as foUows 



9± - ± 



— arctan 

3 \ 9± 



(11) 
(12) 

(13) 



By substituting the mapping transformation (|4]) into the 
relevant expression for the partition function ([3]) one es- 
tabUshes a simple mapping relationship between the par- 
tition function of the spin-1/2 Ising-Hcisenberg model on 
the triangulated Kagome lattice and the partition func- 
tion of the corresponding spin-1/2 Ising model on the 
Kagome lattice 



Z(/3, Jh, a, Ji) = A2^/3Zkag(/3Jf„). 



(14) 



This relation actually completes our exact calculation of 
the partition function as the corresponding exact result 
for the partition function of the spin-1/2 Ising model on 
the Kagome lattice is well known At this stage, 
exact results for other thermodynamic quantities follow 



straightforwardly. The Helmholtz free energy of the spin- 
1/2 Ising- Heisenberg model on the triangulated Kagome 
lattice may be connected to the Helmholtz free energy of 
the corresponding spin-1/2 Ising model on the Kagome 
lattice (.Fkag = — In-Zkag) through the relation 



2iV/3"^ln^/3. 



(15) 



The connection between the internal energy of spin-1/2 
Ising-Heisenberg model on the triangulated Kagome lat- 
tice and that one of the corresponding spin-1/2 Ising 
model on the Kagome lattice can readily be obtained 
by differentiating the logarithm of Eq. with respect 
to the inverse temperature /3 yielding 



dlnZ 



2Nd\nA 




U 



Note that the exact result for the internal energy of the 
spin-1/2 Ising model on the Kagome lattice is well known 
(see for instance Ref. [IH) and thus, our exact calculation 
is essentially completed by introducing the parameters 
Wi and W2 that denote inverse temperature derivatives 
of the functions Vi and V2 given by Eqs. ^ and pH]) 
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The entropy can be now easily calculated either from the 
basic thermodynamic relation between Helmholtz free en- 
ergy and internal energy T = U — TS^ or as a negative 
temperature derivative of the free energy (jlSp . Both pro- 
cedures yield the following closed-form relation for the re- 
duced entropy per one site of the original spin- 1/2 Ising- 
Heisenberg model on the triangulated Kagome lattice 



0,0 



S 



1 

3iV 



■InZ, 



kag 



PW2 / Uk: 



, kag 



(19) 



It is quite obvious from the above formula that the re- 
duced entropy of the spin- 1/2 Ising-Heisenberg model 
on the triangulated Kagome lattice can be expressed 
in terms of the known exact results for the partition 
function^'* and internal energy'^^ of the corresponding 
spin-1/2 Ising model on the Kagome lattice. It is worthy 
to note, moreover, that the zero-field specific heat can 
also be simply obtained as a temperature derivative of 
the internal energy (fTB]) . but the final expression is too 
cumbersome to write it down here explicitly. 



III. RESULTS AND DISCUSSION 

Before proceeding to a discussion of the most interest- 
ing results, it is quite useful to realize that all final results 
derived in the foregoing section are invariant under the 
transformation Ji — > — Jj. For this reason, it is very con- 
venient to set the absolute value of the Ising interaction 
I Ji| as the energy unit and to define two dimensionless 
parameters fceT/l Ji| and Jn/lt/il reducing the number of 
free parameters. The former dimensionless parameter is 
then proportional to a relative size of the temperature, 
while the latter one determines a relative strength of the 
intra-trimer Heisenberg interaction Jh with respect to 
the monomcr-trimcr Ising interaction Jj. 

First, let us take a closer look at the ground-state 
behavior. It is quite evident that the spin-1/2 Ising- 
Heisenberg model on the triangulated Kagome lattice 
exhibits spontaneous long-range ordering if and only if 
the corresponding spin-1/2 Ising model on the simple 
Kagome lattice is spontaneously long-range ordered as 
well. Accordingly, the ground-state phase diagram can 
readily be obtained from a comparison of the effective 
interaction parameter /3J^fg given by Eq. ^ with the 
critical point of the spin-1/2 Ising model on the Kagome 
lattice /3cJkag = ln(3 + 2\/3) [^c = l/(fcBTc), Tc is 
the critical temperature] j22, In the consequence of that, 
the ground state is spontaneously long-range ordered if 
PJ^g > ln(3 -I- 2-^/3), while it becomes disordered spin 

liquid state as long as /3J^^ < ln(3 -I- 2^/3). The phase 
boundary between ordered and disordered ground states, 
which follows from the zero-temperature limit of the 



DISORDERED 
SPIN LIQUID 
PHASE 



fif''<fiJ, =ln(3+2V3)' 

' kag ' c kaf ' ' 




FIG. 2: Ground-state phase diagram in the A — Jh/| Jij plane 
separating the spontaneously ordered and disordered phases. 



mapping parameter ([8]), is shown in Fig. [2] and can be 
expressed through the following analytical condition 



Jh/|Ji| = -2/(2-f A) 



(20) 



It can be easily understood that the disordered spin 
liquid state appears as a result of geometric frustration, 
which comes into play provided that there is a sufficiently 
strong antiferromagnetic Heisenberg interaction. Indeed, 
the ground state is a simple ferromagnetic or ferrimag- 
netic spin arrangement for any Jh/|^i| greater than the 
boundary value (j20p depending on whether Ji > or 
Ji < 0, otherwise it becomes the disordered spin liquid 
state. Apparently, the greater is the anisotropy constant 
A, the weaker antiferromagnetic Heisenberg interaction 
is needed to destroy the spontaneous long-range ordering. 

Now, let us make few remarks on a nature of pos- 
sible ground states. The ferromagnetic as well as fer- 
rimagnetic ordered states can be characterized through 
the same classical spin ordering as previously described 
by the analysis of the spin-1/2 Ising model on the tri- 
angulated Kagome lattice. As far as the disordered 
spin liquid state is concerned, however, there is a fun- 
damental difference between the disordered ground state 
of the spin-1/2 Ising-Heisenberg model with any A 7^ 
and its semi-classical Ising limit A = 0, respectively. 
As a matter of fact, considering the frustrated regime 
>^h/|'/i| < —1 and setting A = into Eq. ([8]) gives the 
zero effective exchange coupling PJ^^g — 0, which means 
that the disordered ground state of the spin-1/2 Ising 
model on the triangulated Kagome lattice is equivalent 
to an ensemble of non-interacting spins, or equivalently, 
to a spin system at infinite temperature. This implies 
that the Ising spins at the monomeric 6-sites are com- 
pletely free to flip and the ground-state degeneracy of 
the disordered spin liquid state can be straightforwardly 
counted by following the argumentation of Loh, Yao, 
and Carlsoni^ The residual entropy per spin is accord- 
ingly So/NrkB = 5 In 72 = 0.4752, since each trimeric 
unit has precisely three different lowest-energy states for 
each possible spin configuration of its three surround- 
ing monomeric spins and the basic unit cell contains two 
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trimeric units and three monomeric spins. 

Contrary to this, the eflFective exchange couphng ^ 
that corresponds to the spin- 1/2 Ising-Heisenberg model 
on the triangulated Kagome lattice with any A 7^ is 
equal to PJ^g = In 2 in the disordered region. Even 
although this value is smaller than the critical value 
/^cAag and the spin system must be consequently dis- 
ordered, its positive and non-zero value indicates a ferro- 
magnetic character of short-range correlations between 
the monomeric Ising spins. It actually turns out that 
the lowest-energy state of each Heisenberg trimer is two- 
fold degenerate if all three surrounding Ising spins are 
aligned alike (all three Ising spins point either 'up' or 
'down'), while there is just one lowest-energy state from 
the ground-state manifold provided that two from three 
Ising spins are aligned alike and the third spin points in 
the opposite direction (remember that the lowest-energy 
state of the trimeric unit is three-fold degenerate irrespec- 
tive of a spin configuration of monomeric spins within the 
Ising model). With regard to this, the residual entropy of 
the disordered spin liquid phase is for the quantum Ising- 
Heisenberg model significantly lower than for its semi- 
classical Ising Hmit So/NrkB = 0.2806 and 0.4752, re- 
spectively, which implies that quantum fluctuations par- 
tially lift a macroscopic degeneracy of the ground-state 
manifold. The most surprising finding stemming from 
our study is that the zero-point entropy of the spin- 1/2 
Ising-Heisenberg model is the same for arbitrary but non- 
zero anisotropy A and this observation suggests that the 
semi-classical Ising limit represents a very special limit- 
ing case of the model under consideration. 

To provide an independent check of the aforementioned 
scenario, we depict in Fig. [3] temperature variations of 
the effective interaction fiJ^^ for several values of the 
ratio Jh/|^i| and two different values of the anisotropy 
constant A = and 1, respectively. It should be men- 
tioned that the particular case with A = 1 qualitatively 
resembles temperature dependences of the mapping pa- 
rameter PJ^^g for the Ising-Heisenberg model with any 
A 7^ 0. The considered model system is spontaneously 
long-range ordered below certain critical temperature if 
and only if the effective interaction parameter /3 J^fg is 

greater than the critical value /3c Jkag = \n{3+2y/3) shown 
in Fig. [3] as a broken line. In agreement with the afore- 
describcd ground-state analysis, the effective interaction 
starts from zero, then gradually increases to some lo- 
cal maximum before it finally goes to zero by increasing 
temperature whenever A = and Jh/I^/iI < —1. On the 
other hand, the effective interaction generally starts from 
the value In 2, then exhibits a temperature- induced local 
maximum before it finally tends to zero whenever A 7^ 
and Jh/|^i| < —2/(2 -I- A). It is also quite evident from 
Fig. [3] that the spin system remains disordered over the 
whole temperature range for any Jh/|>/i| < —2/(2 + A), 
since the effective interaction /3J^ never crosses the 
critical value /3c<^kag that is needed to invoke a spon- 
taneous ordering. This observation would suggest that 
the Ising-Heisenberg model in question cannot exhibit a 
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FIG. 3: Temperature dependences of the effective interaction 
PJkag for several values of the ratio Jh/| Ji| and two different 
anisotropy constants: (a) A = 0; (b) A = 1. Broken line 
shows the critical point /3cJkag ~ ln(3 4- 2\/3) of the corre- 
sponding spin-1/2 Ising model on the Kagome lattice. 



temperature-induced reentrant phase transition from the 
disordered towards the spontaneously ordered phase re- 
gardless of whether A = or A 7^ 0. 

At this point, let us proceed to a discussion of the 
finite-temperature phase diagram. The critical temper- 
ature of the spin-1/2 Ising-Heisenberg model on the tri- 
angulated Kagome lattice can easily be calculated from 
the critical condition /3c^kag = ^^(^ + 2\/3), which en- 
sures that the corresponding spin-1/2 Ising model on the 
Kagome lattice is precisely at the critical point. Note 
that this critical condition is essentially equivalent to a 
graphical solution that finds points of intersection be- 
tween the temperature dependence of the effective inter- 
action PJ^g and the critical point /JcJ'kag — ln(3 + 2\/3) 
of the spin-1/2 Ising model on the Kagome lattice (see 
Fig. [3]). The dependence of dimensionless critical tem- 
perature ksTc/lJil on the ratio J-a/\Ji\ is displayed in 
Fig. |4] for several values of the anisotropy parameter A. 
As one can see, the critical temperature monotonically 
decreases by decreasing the ratio Jh/|Ji| until it van- 
ishes at the ground-state phase boundary ^U\i whenever 
the easy-axis exchange anisotropics A < 1 are assumed. 
By contrast, the interesting non-monotonic dependence 
of the critical temperature may be observed for the easy- 
plane anisotropics A > 1 presumably due to a compe- 
tition between the easy-axis Ising interaction Ji and the 
easy-plane Heisenberg interaction Jh(A) before the criti- 
cal line finally tends to zero-temperature limit consistent 
with the ground-state boundary (|20p . Several limiting 
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FIG. 4: The dimensionless critical temperature k^Tcl\J\\ as 
a function of the ratio Jh/| Ji| for different values of the ex- 
change anisotropy A. Fig. |41d shows a detail of the finite- 
temperature phase diagram, where the critical temperature 
tends towards the zero temperature. 



cases of the model under investigation can be checked. 
First, the critical line of the spin- 1/2 Ising model on 
the triangulated Kagome lattice originally reported by 
Zheng and Sun^l, and later rederived by Loh, Yao, and 
Carlsonp^ is recovered on assumption that A = 0. By 
assuming Jh/|^i| = 0, on the other hand, all critical lines 
meet at a common critical point 



1 



•^il 2 In ( ^3 + 2V3 + ^2 + 2^3 



0.3154, (21) 



which is consistent with the critical temperature of the 
spin- 1/2 Ising model on the decorated Kagome lattice. 
It is also quite interesting to check the asymptotic be- 
havior of the critical temperature achieved in the limit 
JYil\Ji\ — > oo that corresponds to the infinitely strong 
ferromagnetic Heisenberg interaction. If A < 1, the crit- 
ical temperature asymptotically reaches the value 



l^il 



1 



In 



2 + V3 + V6 + 4Vl) 



0.5021, 



(22) 



while for A > 1 it strikingly decreases down to one third 
of this asymptotic value 



1 



-^il 31n(2 + V3 + v/6 + 4V3 



0.1673. (23) 



For the particular case A = 1, the critical temperature of 
the spin-1/2 Ising-Heisenberg model on the triangulated 
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FIG. 5: Thermal variations of the entropy for different val- 
ues of the ratio Jh/|Ji|, which are indicated by the num- 
bers associated with each line, and two different values of the 
anisotropy constant: (a) A = 0; (b) A = 1. The insert shows 
the low-temperature behavior for special cases, where the en- 
tropy tends towards its residual values. 



Kagome lattice acquires in the limit Jh/|^i| oo the 
asymptotic value feB^c/l Ji| = 0.4285. 

Next, let us turn our attention to temperature depen- 
dences of some thermodynamic quantities. Fig. [5] de- 
picts temperature variations of the entropy calculated 
for several values of the ratio Jh/|>/i| and two differ- 
ent values of the anisotropy A = and 1, respectively. 
As it can be clearly seen, the standard S-shaped de- 
pendence with a weak energy-type singularity located 
at critical points of the order-disorder phase transitions, 
which are schematically shown in Fig. O as open circles, 
gradually shifts towards zero temperature upon strength- 
ening the antiferromagnetic Heisenberg interaction un- 
til the weak singularity completely disappears from the 
entropy dependence whenever Jh/|Ji| < —2/(2 + A). 
Under this circumstance, the displayed thermal depen- 
dences tend asymptotically towards the zero-point en- 
tropy So/NTkB = 0.2806 or 0.4752 depending on whether 
A or A = 0, respectively. The low-temperature be- 
havior of the entropy, which is for better clarity shown 
in the insert of Fig. [51 thus provides an independent con- 
firmation of the macroscopic degeneracy that appears 
within the disordered spin liquid ground states. 

Finally, let us conclude our analysis of thermodynamics 
by exploring temperature dependences of the zero-field 
specific heat. For illustration, some typical thermal vari- 
ations of the specific heat are plotted in Figs. [S] and [7] for 
several values of the ratio Jh/|>/i| and two different val- 
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ues of the exchange anisotropy A = and 1, respectively. 
The upper panel in both figures depicts the particular 
cases with the spontaneously long-range ordered ground 
state, whereas the lower panel shows the particular cases 
with the disordered spin liquid ground state. In accor- 
dance with the above statement, all temperature depen- 
dences of the heat capacity shown in both upper panels 
display a logarithmic singularity from the standard Ising 
universality class, which is associated with the continu- 
ous (second-order) phase transition between the sponta- 
neously ordered and disordered phases. Even although 
the specific heat curves displayed in Figs. [H] and [7| have 
several common features, there are nevertheless a few 
important differences. In both cases, the marked round 
maximum is superimposed on the high-temperature tail 
of the specific heat singularity by considering a strong 
ferromagnetic intra-trimer interaction Jh/I^/iI 2> 1 (see 
Figs, and [7^). However, it can be also clearly seen 
that this round maximum is much more pronounced in 
the Ising model with A = than in the Ising-Heisenberg 
model with A = 1. It is worthwhile to remark, more- 
over, that the round maximum in the high-temperature 
tail of the specific heat can be also observed when taking 
into account the antiferromagnetic intra-trimer interac- 
tion Jh/|^i| < (see Figs. [Bbc and [7bc), but in this 
particular case, the robust hump develops only if the ra- 
tio Jh/|>/i| is close enough to the boundary value (PO)) . 
In this parameter space, the most important difference 
between the specific heat curves of the Ising and Ising- 
Heisenberg models lies in an appearance of the additional 
marked shoulder, which appears in the low-temperature 
tail of the heat capacity in the latter model only (Fig.[2b), 
whereas this feature is obviously missing in the relevant 
dependence of the former model (Fig. [51;). 

Last but not least, let us comment on temperature 
dependences of the specific heat by considering some 
typical cases that correspond to a spin system with the 
disordered spin liquid ground state (the lower panels in 
Figs. [Hand [7]). Apparently, the temperature dependence 
with a single maximum is obtained by selecting the ratio 
>/h/|'/i| exactly from the ground-state boundary (j^D)) be- 
tween the spontaneously ordered and disordered phases 
(see Figs. [6]i and[7ti). Note that this maximum is nev- 
ertheless much less symmetric for the Ising-Heisenberg 
model (Figs.[7ti) than for its Ising counterpart (Figs.[njl), 
while both models exhibit a quite similar uprise of the 
additional low-temperature shoulder upon a consecutive 
sfight decrease of the ratio Jh/|=/i|. As evidenced by 
Fig. [7j5, the specific heat of the Ising-Heisenberg model 
gradually looses its irregular shape upon further decrease 
of the interaction parameter Jh/|^i|. It is noteworthy, 
moreover, that the specific heat curve with two marked 
rounded maxima appear by assuming a sufficiently strong 
antiferromagnetic intra-trimer interaction Jh/|^i| ^ — 1 
(Figs.[6f and[7f), whereas the more negative is the ratio 
Jh/|Ji|, the more evident is the double-peak structure 
of the zero-field specific heat curve. However, it should 
be also pointed out that the temperature dependence 



with the double-peak specific heat emerges in the Ising- 
Heisenberg model with A = 1 for much smaller geometric 
frustration (i.e. less negative ratio Jh/|^i|) than for the 
Ising model with A = 0. Notwithstanding this obser- 
vation, the specific heat curves with two well-separated 
maxima might be expected for the series of polymeric co- 
ordination compounds Cu9X2(cpa)6 for which a rough es- 
timate of the ratio \Jaa/ Jab\ ~ 40 has been made accord- 
ing to the experimentally measured susceptibility datai^ 



IV. CONCLUDING REMARKS 

The present article provides a survey of exact analyt- 
ical results for the spin- 1/2 Ising-Heisenberg model on 
the triangulated Kagome (triangles-in-triangles) lattice, 
which has been proposed in order to shed light on a frus- 
trated magnetism of the series of Cu9X2(cpa)e polymeric 
coordination compounds. The model under investiga- 
tion has exactly been solved with the help of general- 
ized star-triangle transformation that establishes a pre- 
cise mapping relationship to the corresponding spin- 1/2 
Ising model on the Kagome lattice. Exact results for 
the ground-state and finite-temperature phase diagrams, 
as well as, several thermodynamic quantities such as the 
Helmholtz free energy, internal energy, entropy, and spe- 
cific heat, have been derived and discussed in detail. 

Even though a theoretical description based on the hy- 
brid Ising-Heisenberg model may not be fully realistic 
for true magnetic materials from the family of copper- 
based coordination compounds (Heisenberg model would 
be of course more adequate), it is quite reasonable to ex- 
pect that the presented exact results illustrate many im- 
portant features of the series of three isostructural com- 
pounds Cu9X2(cpa)6. It is really plausible to argue that 
the suggested model correctly takes into account quan- 
tum effects closely connected with the stronger intra- 
trimer Heisenberg interaction Jh , whereas the Ising char- 
acter of the weaker monomer-trimer interaction Ji might 
be regarded as at least rather reasonable first-order ap- 
proximation. This approximation should be altogether 
acceptable especially in the highly frustrated region 
J'a/\Ji\ —2/(2 -|- A), where the 'Ising' spins located at 
the monomeric sites are completely free to flip without 
any energy cost owing to a strong geometric frustration 
produced by the trimeric 'Heisenberg' spins. From this 
perspective, it is quite tempting to conjecture that the 
zero-field specific heat of CugX2(cpa)6 compounds should 
exhibit a notable temperature dependence with two out- 
standing round maxima, which are well separated from 
each other because the antiferromagnetic intra-trimer in- 
teraction is roughly two orders of magnitude stronger 
than the ferromagnetic monomer-trimer interaction. 

In agreement with our expectations, it also has 
been demonstrated that quantum fiuctuations intro- 
duced through the non-zero exchange anisotropy A help 
to destroy spontaneous (ferromagnetic or ferrimagnetic) 
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long-range ordering. As a matter of fact, it can be read- 
ily understood from Eq. ([20|) that the increase in the 
anisotropy parameter A suppresses a strength of the an- 
tiferromagnetic intra-trimer interaction that is needed to 
prevent spontaneous ordering. However, the most sur- 
prising finding to emerge from the present study closely 
relates to a substantial decline of the residual entropy 



of disordered spin liquid state, which appears on as- 
sumption that there is arbitrary but non-zero exchange 
anisotropy A. It actually turns out that the quantum 
fluctuations partially lift a rather high macroscopic de- 
generacy of the disordered spin liquid state and conse- 
quently, the zero-point entropy of the quantum Ising- 
Heisenberg model So/NTk^ = 0.2806 is for any A 7^ 
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significantly lower than the zero-point entropy of the 
semi-classical Ising model So/Nj^kB — 0.4752 achieved 
in the limiting case A = 0. From this point of view, 
another interesting question arises whether or not the 
residual entropy of the spin- 1/2 Heisenberg model on the 
triangulated Kagome lattice will be completely removed 
by the order-from-disorder mechanism^^ when account- 
ing for the Heisenberg character of the monomer-trimer 
interaction as well. In order to clarify this unresolved 
issue, it is necessary to perform a rather extensive exact 
numerical diagonalization of the full Heisenberg Hamilto- 
nian, or, alternatively, it would be also possible to use the 
exact solution of the Ising-Heisenberg model as a starting 
point of more general perturbative treatment. 

Note added 

Shortly after our manuscript has been submitted for 
publication, Yao et al^ has completed a similar work 
dealing with the spin-1/2 Ising-Heisenberg model on the 
triangulated Kagome lattice. The main difference be- 
tween our procedure and the one of Yao et al.^'^ lies in 
the way of calculating the energy spectrum of Heisenberg 
trimer and its three enclosing Ising spins. In our proce- 
dure, we have firstly performed an exact diagonalization 
of the Heisenberg trimer in some generally non-uniform 
local field produced by three surrounding Ising spins and 
then, we have considered two symmetry-distinct config- 
uration of the enclosing Ising spins in order to establish 
a precise mapping relationship with the corresponding 



spin-1/2 Ising model. In the procedure developed by Yao 
et al.,^'' the authors first consider two symmetry-distinct 
configurations of the enclosing Ising spins and then, they 
diagonalize the Heisenberg trimer for both symmetry- 
distinct configurations of the surrounding Ising spins. In 
a such way, they were able to avoid rather cumbersome 
expressions for the roots of cubic equations, which other- 
wise occur when diagonalizing the Heisenberg trimer in 
some generally non-uniform local magnetic field. Despite 
this difference, both procedures turn out to be completely 
equivalent as it could be easily checked from a comparison 
of the effective mapping parameter (iJ^^ calculated from 
our Eqs. (8-13) and Eqs. (7-8,13-14) of Yao et al.^ re- 
spectively. Exact results presented in both manuscripts 
for the ground-state and finite-temperature phase dia- 
grams, Helmholtz free energy, internal energy and en- 
tropy are consequently in a very good accordance. 
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